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We performed first-principle calculations based on the supercell and cluster approaches to inves- 
tigate the neutral Al impurity in smoky quartz. Electron paramagnetic resonance measurements 
suggest that the oxygens around the Al center undergo a polaronic distortion which localizes the 
hole being on one of the four oxygen atoms. We find that the screened exchange hybrid functional 
successfully describes this localization and improves on standard local density approaches or on hy- 
brid functionals that do not include enough exact exchange such as B3LYP. We find a defect level at 
about 2.5 eV above the valence band maximum, corresponding to a localized hole in a O 2p orbital. 
The calculated values of the g tensor and the hyperfine splittings are in excellent agreement with 
experiment. 



I. INTRODUCTION 



Silicon dioxide is among the most commonly encoun- 
tered substances in both daily life and in electronic ap- 
plications. It is highly abundant in nature in the form of 
crystalline quartz and can be grown extremely pure by 
experimental techniques. A prominent defect is the neu- 
tral [A104]° center, which has been identified by its elec- 
tron paramagnetic resonance (EPR) signatur#HSl Here, 
an Al'^'*'-ion substitutes a Si-ion, leaving an unpaired elec- 
tron at one of the four oxygen atoms adjacent to the 
Al center. The corresponding localized spin gives rise 
to an ESR signal. This defect is believed to cause the 
smoky coloration of quartz crystals^ ^. On theoretical 
level, this 'classical' model has been confirmed by clus- 
ter calculations on unrestricted Hartree-Fock or hybrid 
functional levePH^i' , yielding a polaronic hole localiza- 
tion and hyperfine coupling constants with ^^Al, ^''O 
and ^^Si in good agreement with experiments. It was 
shown that the defect undergoes local symmetry break- 
ing and Jahn- Teller reconstruction with one oxygen atom 
relaxing away from the Al ion. At the other hand, den- 
sity functional theory (DFT) calculations predicted that 
the hole and spin are delocalized over all four oxygen 
ato ms an d no symmetry breaking between the oxygen 
site^i^IIll Similar predictions occured for other oxide 
wide-gap semiconductors, where DFT on the level of the 
local density (LDA) or generalized gradient approxima- 
tion (GGA) level fails to correctly predict the polaronic 
hole localization on oxygen, such as in case of cation va- 
cancies in ZnOC^HISl^ mOp^, TiO^ and acceptor im- 
purities in GalsP, InaOs or SnaO;^. The observed 
delocalization of the hole arises from the residual one- 
electron self-interaction from the Hartree energy, which is 
insufficiently canceled by the exchange-correlation term 
of LDA and GGA type functionals, thus promoting arti- 
ficial delocalizatioEpil. 

Similarly, the self-interaction also contributes to the un- 
derestimation of experimental band gaps by predicting 
valence (conduction)bands at too high (low) energies. 
Indeed. Mauri et aZp2! showed that self-interaction cor- 



rected LDA^'^ calculations, where the contributions due 
to self-interaction are explicitly substracted in the energy 
expression, favor a distorted geometry and trapped hole. 
On the other hand, Hartree-Fock exchange completely 
cancels self-interaction contributions and inclusion of HF 
exchange in the exchange-correlation functional typically 
can compensate for the shortcomings of LDA and GGA. 
However, pure HF results in a gross over-estimate of the 
band gap. It is interesting that some of these hybrid 
functionals, such as B3 LYP, do not give the full hole lo- 
calization in SiOdiSlMI, This suggests that a minimum 
amount of exact exchange is needed to give this result 
correctljJ^H. 

In this paper, we investigate the [AIO4]*' center using 
periodic boundary conditions and the screened-exchange 
(sX-LDA) hybrid functional, which includes a screened 
version of the full exact exchange. It is essentially self- 
interaction free for all electron spacings less than the 
screening length and improves on the predicted band 
gapg2^H2Z!, We show that sX-LDA can restore the local- 
ization of the polaronic hole on one oxygen atom. We fur- 
ther provide calculated EPR parameters, which we find 
to be in good agreement with experiment. To the best of 
our knowledge, calculated Lande g tensors of the [AIO4]'' 
center have not been reported so far. 



II. METHOD 

The calculations were performed in the frame of spin- 
polarized density functi onal t heory (DFT) using the hy- 
brid functional sX-LDA'^SlSfil^ which has been recentlji^H 
implemented in the planewave code CASTEP^^. Here, 
the self-energy of an electron in the crystal is approxi- 
mated by a combination of a short-range Hartree-Fock 
exchange-type term 
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and a long range, local density dependent. LDA-type 
term 

This range-separation into short and long range terms 
is similar to that in other func tional s like the Heyd- 
Scuseria-Ernzerhof (HSE) versioiP^'^. In theory, the 
Thomas- Fermi wave vector kj-p, which depends on the 
average charge density, is used as the inverse screening 
length ks- In practice, the TF wavevector needs to be 
chosen carefully in terms of the density of valence elec- 
trons. Here we use a fixed value of 0.76 bohr, which works 
well for sp semiconductors. 

The calculations were done in two steps: In the first step, 
we used periodic boundary conditions to model the defect 
in the solid by a 2x2x2 supercell of the common (hexago- 
nal) 9 atom unit cell of a-Si02 with an Al atoms substi- 
tuting one of the 24 Si atoms, and optimized the atomic 
positions while keeping the lattice constants at the ex- 
perimental values. The integrals in reciprocal space were 
approximated by the values at the Baldereschi k-point for 
hexagonal latticeJ^and we used OPIUM'^'' pseudopoten- 
tials with a cutoff energy of 800 eV to model the valence 
electrons of Si and O. 

In the second step, we calculated the electron param- 
agnetic resonance (EPR) parameters with the Quantum 
Chemistry code ORCA54, which particularly aims at the 
spectroscopic properties of open-shell molecules. We 
used a 33 atoms cluster, which we obtained by cutting 
out the first three atomic shells surrounding the Al cen- 
ter from the previously relaxed 72 atom supercell and 
passivating the dangling bonds of the outmost 12 oxy- 
gen atoms by hydrogen. The atoms were represented by 
all-electron def2-TZVP sets from the Karlsruhe groups. 
The exception are the oxygen atoms directly adjacent to 
the Al-center that were described by Barone's EPR-IlPI 
basis set, which was specifically designed for the calcula- 
tion of EPR properties. As ORG A does not offer screened 
hybrid functionals, we chose to use to use Becke's 'Half 
and Half (BHandHLYP)^^ hybrid functional 

instead. 
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FIG. 1: (Golor online) Partial density of states of a 72 
atom supercell of Si02:Al from sX-LDA calculations. 
The dashed lines represent the energy of the valence 
band maximum and the conduction band minimum. 




FIG. 2: (Golor online) Gorresponding plot of the spin 
density of Si02:Al around the Al center. The system 
undergoes a polaronic distortion and the spin is 
localized almost exclusively at one of the four oxygen 
atoms adjacent to the Al impurity. Some atoms are 
hidden by the spin density isoplane. 



III. RESULTS AND DISCUSSION 

Fig. [T] shows the calculated partial density of states 
(PDOS) for a 72 atom supercell of Q:-Si02 containing 
a single Al center. Glearly, the impurity introduces 
a single spin-polarized defect state at an energy of 
2.5 eV above the valence band maximum. The defect 
state is predominantly of O 2p character with small 
contributions from O 2s and Si 3p mixed in. We find 
that the overall magnetic moment of 0.99 hb of the 
system is localized almost exclusively at one single 



oxygen atom, refer to the spin density plot in Fig. [2] 
Together with the localized hole, our sX-LDA approach 
yields an asymmetric distortion of the geometry sur- 
rounding the Al atom. We find an Al-0 bond length 
of 1.92 A for the oxygen atom carrying the localized 
hole, corresponding to an bond elongation of 13% 
compared to the other three Al-0 bonds (1.69-1.7lA). 
This prediction is in reasonable agreement with esti- 
mation from EPR measurement^l2% bond elongation). 

We note that the degree of spin localization of depends 
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TABLE I: Calculated principle values of the g-tensor and the hyperfine matrix of the [A104]'^ center from 
calculations on a Al(04SiH3]4 cluster, and comparisons with previous theoretical from other authors and 
experiments. Following Ref.El the oxygen labels O* and O^'^''^ refer to the oxygen atom with the hole and the three 
remaining oxygen atoms adjacent to the aluminium center, respectively. The values for ^^Si are given for the Si 
nucleus with he largest contribution. 



This work 



BBl 



Exr 



Spin 
(in ^b) 



O* 

01.2.3 



0.93 
0.01-0.05 



1.04 
<0.01 



0.81 
0.00-0.20 



g tensor 



gi 

g2 
g3 

Siso 



2.0031 
2.0093 
2.0412 
2.0178699 



2.0024±0.0003 
2.0091±0.0003 
2.0614±0.0003 



170* 

hyperfine 
coupling 
matrix 
parameters 
(in G) 



Ai 
As 
As 

A. 

■^JSO 

Bi 

B2 
Bg 



-119.225 
22.544 
22.86 
-24.607 
-94.618 
47.151 
47.467 



-128.6 
11.6 
13.6 
-34.5 
-94.1 
46.1 
48.1 



-109.46 
23.48 
23.76 
-20.74 
-88.72 
44.22 
44.50 



-111.0 
15.2 
17.8 
-26.0 
-85.0 
41.2 
43.8 



2^A1 
hyperfine 
coupling 
matrix 
parameters 
(in G) 



Ai 
Aa 
As 

A. 

Bi 

B2 
Bg 



-6.236 
-5.972 
-4.953 
-5.721 
-0.516 
-0.251 
0.767 



-5.1 
-5.1 
-4.2 
-4.8 
-0.3 
-0.3 
0.6 



-7.09 
-6.99 
-5.90 
-6.66 
-0.43 
-0.33 
0.76 



-6.2 
-6.1 
-5.1 
-5.8 
-0.4 
-0.3 
0.7 



hyperfine 
coupling 
matrix 
parameters 
(in G) 



Ai 

A2 

As 

A. 

Bi 

B2 

Bs 



8.744 
9.347 
9.524 
9.205 
-0.403 
0.142 
0.319 



15.5 
17.4 
17.8 
16.9 
-1.4 
0.5 
0.9 



10.11 
10.57 
10.77 
10.48 
-0.37 
0.09 
0.29 



10.8 
11.4 
11.6 
11.3 
-0.5 
0.1 
0.3 



^ Obtained with an Al(OSiH3)4 cluster. 
^ Obtained with an AlOieSiisHae cluster. 



on the included portion of Hartree-Fock exchange in 
the functional. Pacchioni et al^ reported that the 
popular B3LYP functional (20% HF exchange) predicts 
a partial spin localization with spin populations of 0.29 
on two oxygen atoms and 0.21 on the remaining two. 
Correspondingly, the Al-0 bond lengths in this case are 
very similar. On the other hand. To et al'J^ found a 
bond elongation of 12% using the BE IK functional (42% 
HF exchange), whereas Pacchioni et al'^veported a 14% 



elongation from their 100% Hartree-Fock calculations. 
Our cluster calculations using the BHandHLYP hybrid 
functional (50% exact exchange) yield spin localization 
and defect geometry very similar to our sX-LDA results. 



In the next step, we calculated the EPR parameters 
for the defect. The EPR spectrum can be modeled by an 
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Hamiltonian 



iIe// = fs.g.B + ^S.A.I, 



(1) 



where a is the fine structure constant. The first term of 
Eq. [l] describes the coupHng of the spin momenta S of 
the unpaired electrons with an external magnetic field B 
by the g tensor 




The second term of Eq. [T] represents the hyperfine cou- 
pling of the electron spin with the nuclear spins can be 
described in terms of the hyperfine matrix A, which can 
be separated into its isotropic and its anisotropic part B 
(related to dipolar interaction) 
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A problem arises from the fact that the g tensor is not 
gauge invariant, so that the results technically depend on 
the choice of origin in the calculation. While the origin 
dependence is usually small as long as sufficiently large 
basis sets are used, we have checked the sensitivity of the 
results with respect to the origin. We found that on GGA 
level the center of the electron charge gives practically 
the same values as calculations using the fully invariant 
IGLO (individual gauges for localized orbitals)'^^ proce- 
dure (which, however, is not rigorous for hybrid func- 
tionals) and used this point as the origin for our cal- 
culations. Table |T] shows the obtained principle values 
of the g-tensor and the hyperfine matrix in comparison 
with previous theoretical results and experiment. In case 
of the g-tensor, our computed values for the two smaller 
principle values 5i=2.0031, (72=2.0093 are in very good 
agreement with the experimental estimations gi „=2.0024 
and (72,0=2.0091. The deviations for the third principal 
value are larger. We find (73=2.0412, which is consider- 
ably smaller than (73 ^=2.0614 as reported by Schnadt et 

For the hyperfine coupling, we find a good prediction 
of our approach for the investigated ^''O, ^^Al and '^^Si 
nucleii. This is particularly true for ^^Al, where our re- 
sults, being the golden mean of the reported values of 
UHpii^ and BBII^^ calculations, are remarkably close 
to experiment. For ■^^Si, our calculations slightly under- 
estimate the experimental values and yield similar results 



to those from To et al. This is entirely due to the un- 
derestimation of the isotropic component of A, A^^^ by 
2G, while the (small) dipolar interaction is reproduced 
accurately. The dominant contribution to hyperfine cou- 
pling is to be expected from the oxygen atom with the 
localized hole. Here, the electron-nucleus interaction con- 
tains a strong dipolar component and the experimen- 
tally obtained hyperfine matrix shows a strong anisotropy 
for the three principal axes. The UHF calculations by 
Pacchioni et al. underestimated the experimental val- 
ues Ai,e=-111.00G, A2,el=15.2G and A3,e=17.8G due 
to a too large isotropic contribution. In contrast, the 
BBIK calculations predicted a considerably weaker Ai^^, 
resulting in a slight underestimation of Ai and overes- 
timation of A2 and ^3 and the opposite behavior for 
the anisotropic matrix. We find the hyperfine parame- 
ters Ai=-119.22G, A2=22.54G and A3=22.86G, which 
overestimate all three of the experimental values. The 
reason is a too strong anisotropic component in our hy- 
perfine matrix, as Tab |T] shows. At the other hand, our 
predicted A^^^ is favorably close to the experimental value 
and a better prediction than the values from the other 
two methods. In fact, the differences between our param- 
eters and those from Pacchioni et al. arise from the differ- 
ent isotropic contributions, while the dipolar interaction 
is close to identical in both cases. The dipolar contribu- 
tion depends on the angular momentum of the unpaired 
electron and hence is roughly proportional to the popu- 
lation of the O 2p orbital and the localization of the spin 
density. To et al have reported that increasing the size 
of the cluster in their case led to a stronger spin local- 
ization and a larger dipolar interaction. We could thus 
argue that our calculations slightly overestimate the lo- 
calization of the electron and the correct spin population 
of the oxygen atom is about 0.9 ^b. This is in line with 
the observation that overall the slightly smaller portion 
of Hartree-Fock exchange in the BBIK functional com- 
pared to BHandHLYP seems to benefit the prediction of 
hyperfine parameters for all investigated nucleii. 



IV. CONCLUSION 

We showed that sufficient inclusion of Hartree-Fock in 
hybrid functionals does improve on DFT supercell cal- 
culations and can correctly describe the polaronic hole 
and the corresponding symmetry distortion of the neu- 
tral Al impurity in a-qaaxtz. The observed localization 
is inherently connected to the reduced self-interaction in 
hybrid functionals, which shows in by its dependence on 
the ratio of HF:DFT exchange. The validity presented 
approach was further shown by a subsequent calculation 
of the g-tensor and the hyperfine coupling matrix using 
a cluster approximation for the sX-LDA geometry, which 
was in good agreement with the conclusions from previ- 
ous EPR measurements. 
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